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Abstract 

We complete the determination of the universal amplitude ratios of two-dimensional perco- 
lation within the two-kink approximation of the form factor approach. For the cluster size 
ratio, which has for a long time been elusive both theoretically and numerically, we obtain 
the value 160.2, in good agreement with the lattice estimate 162.5 ± 2 of Jensen and Ziff. 



1. Universal combinations of critical amplitudes represent the canonical way of encoding the 
universal information about the approach to criticality in statistical mechanics [T]. While critical 
exponents can be determined working at criticality, amplitude ratios characterize the scaling 
region around the critical point. They carry independent information about the universality 
class and their determination is in general theoretically more demanding. Field theory is the 
natural framework in which to address the problem, but the usual perturbative approach is not 
helpful if one has to work far below the upper critical dimension d c . 

For the best-known example of a geometric phase transition, namely isotropic percolation 
{d c = 6) [2], it was shown in [3] how the field theoretical computation of universal amplitude 
ratios in two dimensions can be addressed non-perturbatively exploiting the fact that percolation 
can be seen as the q — > 1 limit of the g-state Potts model [3], and that the latter is integrable 
even away from criticality, in the scaling limit for q < 4 [5]. Starting from the exact 5-matrix 
[5] one can compute the Potts correlation functions, and from them the amplitude ratios, using 
the form factor approach [3J. 

This programme was completed in [3] for q = 2, 3, 4, recovering the known Ising results and 
obtaining new predictions for the three- and four-state Potts model. For percolation, however, 
only partial results were obtained, because the determination of some amplitudes above the 
percolation threshold p c involves the solution of a functional equation which in [3] could not be 
solved for generic values of q, in particular for q = 1. In this situation, it was observed in [3] 
that a simple parabolic extrapolation to q = 1 of the results obtained at q = 2, 3, 4 produced for 
the percolation ratios results compatible with the available numerical estimates. In particular, 
for the ratio of cluster size amplitudes below and above p c the extrapolated value 74 essentially 
coincided with the central value of the most recent estimate then available [6]. 

On the other hand, the status of the numerical results for the cluster size ratio (associated 
to the Potts susceptibility ratio) was at the time particularly controversial, different authors 
having obtained over the years values which spanned two orders of magnitude [T]. Following 
the appearance of [3], Jensen and Ziff communicated a new, very accurate lattice determination 
of this ratio, essentially coinciding with the value 162.5 ± 2 finally published in [7]. A possible 
explanation for such a large discrepancy, other than the failure of the extrapolation (which 
appeared to work for other ratios), was discussed in [5], but was ruled out by the full numerical 



confirmation the prediction of [31 [8] for the universal ratios at q = 3 received_| in [9j [TO] (see also 



In this paper we provide the only piece of analytic information missing in [3], namely the 
solution of the functional equation at q = 1, and show that it leads to results for the percolative 
universal ratios in complete agreement with the most recent lattice estimates. In particular, this 
confirms that the only problem with the extrapolated value for the cluster size ratio was the 
extrapolation itself. When comparing field theoretical and lattice results (Table 2 below) it must 
be taken into account that, with few exceptions, the former are themselves not exact, since they 

1 For the case q = 4, which is plagued by logarithmic corrections to scaling [121 113] , the issue of the precise 
comparison between field theoretical and lattice results for the universal ratios appears still open [141 [S] [TS] . 
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are obtained truncating the spectral series for correlation functions to the two-kink contribution. 
The remarkable accuracy of this two-particle approximation is however well known (see e.g. the 
comparison with the Ising exact results in [3]), and is further illustrated by this case. 

2. In [3] the determination in the two-kink approximation of the low-temperature spin-spin 
correlation function of the scaling g-state Potts model (q < 4) was reduced to that of a function 
£lq(0) entering the two- kink form factor of the spin field. This function is characterized by the 
following properties [3]: 

i) is a meromorphic function of 9 whose only singularity in the strip Im 9 6 (0, 2ir) is a simple 
pole at 9 = in with residue 

Res e=l M q (0) = i^—M, (1) 
q — 1 

where M denotes the Potts spontaneous magnetization; 

ii) is solution of the functional equations 

n q (9) = n q (-e), (2) 

2 cos — sinh X9 tl q (9) = sinh A(nr + 9) Q q (2i-K -9)- sinh \{m - 9) tl q (2iir + 9) , (3) 
3 

with the asymptotic behavior 

" '2 



fl q (9) ~ exp 



A-l 
3 



+oo, (4) 



where A parameterizes q according to the relation yjq = 2sin(-7rA/3). 

For q < 3, where the Potts scattering theory possesses no bound states, the properties i) and 
ii) uniquely identify £l q (9), and then the spin field of the scaling Potts modej§. 

In [3] Q q (9) was determined only for q = 2,3,4, where it takes a simple form. We now show 
that £li(9) can be obtained taking a mathematical detour in the sine-Gordon model. For the 
latter, which is defined by the action 

Asg = J d 2 x ^ {d u tp) 2 + ii cos pip, (5) 

Lukyanov computed in [18] the soliton-antisoliton form factors 

KeM = <0|e^°) 1^(01)^(02)), St = ±1, 8,82 = -1, (6) 

obtaining a result that in our notationg§| reads 

F* { e) = -(e ia ^)^p,A%{9), (7) 

A a ± {9) = e T ^- e \e^™I a {-9) + J a (0 - 2m)}, (8) 



2 See [161 II 7J for the correspondence between fields and solutions of the form factor equations in integrable field 
theory. 

3 In particular, switching from Lukyanov's notations to ours involves the replacements 8 — > —9, £ — > !;/tt, 
a — > /3a. 
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where £ = 7r/3 2 /(87r — /3 2 ), Fq(9) is a function on which we comment below, and I a (9) is specified 
for real values of 9 and a € (— | — |, i) as 



ho ° da; 



-x -\ h ivr I e V ? 

2 



-2a he 



with 



and 



cosh ( 



exp 



oo ^ sinh ^ ( 1 



t sinh i sinh P- 

Z7T 



— (in — 9) 



2n 



C = - exp 

4 



00 dt 

T 



sinh 2 4 sinh 4(1-- 
sinh t sinh 



(9) 



(10) 



(11) 



These results were presented in [18] within a framework known as free field representation, 
which differs from the usual approach to form factors based on functional relations. Of course, 
this latter approach can be adopted also for the matrix elements ([7]), using as input the sine- 
Gordon S"-matrix and the fact that the soliton is semi-local with respect to e^^, with a semi- 
locality phase e 2 "™. The corresponding functional equations then read [19] 



KeM = S T (9)F £ a 2£l (-9) + S R (6)F? ie2 (- 



where 



F £ a l£2 (9 + 2in) = e 

St(0) = 
Sr(9) = 



sinh^ 



sinh j(0 — in] 
sinh ^ 



-0), 



sinh I ( 



ITT) 



(12) 
(13) 



(14) 
(15) 



are the transmission and reflection amplitudes; the explicit form in the present notations of S(9) 
and Fq(9) can be found for example in [20], but here we only need to know that 



F (9) = S(9)F (-9), 



F (9 + 2in) = F (- 



(16) 



This implies in particular that f)13[) is automatically satisfied by ([7]). Since A± are meromor- 
phic functions of 6, also I a , as a linear combination of and A°l with entire coefficients, is 
meromorphic. It particular, analyticity implies that the property 



i a (e) = i a (-e), (17) 

which for real 9 is apparent in ([9]), extends to the whole complex plane. Requiring (|12h leads 
now to the equation 

2cos + 2vra^ sinh y I a (9) = sinh j(in-rj9) I a (2m- 9) -sinh j(in + rj9) I a (2in + 9) , (18) 
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with rj = 1. Making the identifications 



Z = j, o = ~fl±|) 



we rewrite (fTBj) as 



7rA 

2 cos — sinh A<9 4(0) = sinh \(iir - r]0) I a (2iTr - 0) - sinh A(ivr + t}0) I a (2i7r + 0) , (20) 
3 

always with rj = 1. On the other hand, this equation coincides with ([3]) when rj = — 1. At g = 1 
(i.e. A = 1/2), however, the sign of r\ becomes immaterial and (|20p exactly coincides with the 
equation satisfied by 

The functional equation (|20p has infinitely many solutions (a solution multiplied by a 2i7r- 
periodic function of 8 is a new solution) and it remains to be seen whether Q with the identi- 
fications ()19p and A = 1/2 can yield the function Oi relevant for the percolation problem. 

From the known asymptotic behavior (see e.g. |19j ) of the form factors (|7J) one deduces that 
I a (6) behaves as exp [(a — ^) 0] as 6 — > +oo, a result which is not obvious from Q but can be 
checked numerically. Comparing with we see that I a behaves asymptotically as f?i provided 
we take a = —1/6, corresponding to the lower sign and k = in (|19p with A = 1/2. 

The value £ = 27r (i.e. A = 1/2) falls in the repulsive regime of the sine-Gordon model in 
which the only singularity of the form factors (|7J) within the strip Im € (0, 2ir) is the annihilation 
pole at 8 = in. Since Fq(6) is free of poles in the strip, the annihilation pole must be carried 
by A±, and then by I a . Any other pole of I a in the strip could not cancel simultaneously in 
and A°L, and then I a (6)\^ =2n possesses a single pole at 6 = in in the strip ImO 6 (0, 2ir), exactly 
as it is the case for Q q (0) in the Potts model. 

Summarizing, the functions I^i/e(0)\^_ 2iT and Qi(0) satisfy the same functional relations, 
have the same asymptotic behavior and the same singularity structure; then we conclude that 
they coincide up to a normalization. Since we know that [19] 



Res e=ln F«_(8) = i(l - e- 2i ™)<e^>, (21) 

we read from ((7J), © and (fT7|) that I a (0) has residue i on the pole. Knowing also that the 
percolative order parameter P (probability that a site belongs to an infinite cluster) is related 
to the Potts magnetization aa3 

P = lim — - — M, (22) 
and recalling ([1]), we conclude that fii(0) has residue iP on the pole, and therefore 

n 1 (e) = p/_ 1/6 (e)| f=ar . (23) 

The values £ = 2-/r, a = —1/6 fall in the range where (|9|) can be used to compute f2i(0) for real 
values of 8, which is sufficient for our purposes. 



4 The relation (|22[) is written incorrectly in [3], see [8]. 
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3. Near the percolation threshold the relations 

S-r+lp-pd-T, (24) 

i-f^p-Pcf, (25) 

P~B(p-p c f, (26) 

{Nc) - A ± |p-p c | 2 -«. (27) 



N 

define the critical amplitudes for the mean cluster size, the correlation length, the order param- 
eter and the mean cluster number per site, respectively; the superscripts ± refer to p c being 
approached from below or from aboveH. Below we consider both the second moment correlation 
length 

2 _ 1 fd 2 x \x\ 2 g c (x) 

?2nd — T Tl2 T~\ — ' V Z8 J 



and the true correlation length £j defined by 

g c (x) ~ e~ |x|/€f , \x\-Kx, (29) 

where g c (x) is the probability that x and the origin belong to the same finite cluster. It was 
shown in [3] that, in terms of the Potts kink mass m, £j is 1/m at p < p c and l/2m at p > p c , 
and that 

A ± = 7 J— -. (30) 

2 V / 3(/ t + ) 2 1 ' 

Defining the amplitude combinations 

= [a(l - a)(2 - a)A+] 1/2 /+, C7 = 4^^^, (31) 

which are universal due to the scaling relations 2 — a = 2u and 2v = 2(3 + 7, (|30p together with 
a = —2/3 imply in particular 

1 1/2 



5* 



40 
27^/3 



0.9248.., (32) 



a result recovered from a lattice computation in [21]. The result for i?^ in the two-kink 
approximation was computed in [3] and compares quite well with the lattice estimate obtained 
from the combination of the data collected in Table 1. 

The result (|23p allows us to complete the two-kink computation of the universal ratios 
involving the amplitudes f^nfr ^- -> All the other necessary information was already given in 
[3] and here we only recall how the results for percolation follow from those for the Potts model 
when q — > 1. 

Consider as an example the cluster size S. This is the limit of the Potts susceptibility divided 
by q — 1 , and the susceptibility is in turn the integral on the plane of the connected Potts spin- 
spin correlator. At T < T c the leading large distance contribution to this correlator is produced 



5 We keep the notation of [3] where ± referred to the high/low-temperature Potts phases; we drop instead the 
tilde used there on percolation amplitudes. 
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Triangular 


Square 


A+ 


-4.37 a 




r+ 


0.0720 6 


0.102 6 


B 


0.780 b 


0.910 6 


^ftnd 


0.520 6 


0.520 6 



Table 1: Lattice estimates of critical amplitudes for site percolation on triangular and square 
lattice. The superscripts a, b indicate Refs. [22j[23], respectively. 



by a two-kink state and is multiplied by q — 1 (the number of two-kink intermediate states in the 
low-temperature spectral sum). There are no other factors of q — 1 since in this phase the spin 
two-kink form factor F[(9) is the product of £l q (0) times another function which is also finite at 
q = 1 (see [3] ) . At T > T c the spin-spin correlator coincides by duality with the low-temperature 
disorder-disorder correlator. The latter receives the leading contribution from a single one-kink 
state weighted by the squared disorder form factor \F^\ 2 = M\F[ (oo)| [3]. As a consequence, 
due to (|22p , also the high-temperature Potts correlator vanishes as q — 1 in the percolation limit 
(the two- kink contribution behaves in the same way). Summarizing, the factors of q — 1 can be 
explicitly isolated and cancel in the computation of the percolative critical amplitudes for the 
cluster size. The same can be shown for the other amplitudes. 





Field Theory 


Lattice 


A+/A- 


1 


l a 


ft/ft 


2 




find/ ft 


1.001 




find/ f 2nd 


3.73 


4.0±0.5 C 


r+/r~ 


160.2 


162.5 ± 2 d 


u 


2.22 


2.23±0.10 e 


Rj 


0.926 


» 0.93 a+b 



Table 2: Universal amplitude ratios in two-dimensional percolation. The field theory results 
in the first two lines are exact, the others are obtained in the two-kink approximation. The 
superscripts a, b, c, d, e indicate Refs. [22j [231 El El El] , respectively. 

The field theoretical results for the complete list of independent ratios involving the am- 
plitudes (I24H27P are summarized in Table [2] together with the most accurate lattice estimates. 
As remarked above, the comparison confirms in particular the effectiveness of the two-particle 
approximated form factor results in integrable field theory. 

6 In [3] the ratio R c — 4(i?^" d ) 2 /U was considered instead of U. The result R c = 1.56 we obtain should be 
compared with « 1.53 following from Table [1] 
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